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Abstract 

Solutions of hydrodynamical equations are presented for an equation of state al- 
lowing for a first-order phase transition. The numerical analysis is supplemented by 
analytical treatment provided the system is close to the critical point. The processes 
of growth and dissolution of seeds of various sizes and shapes in meta-stable phases 
(like super-cooled vapor and super-heated liquid) are studied, as well as the dynam- 
ics of unstable modes in the spinodal region. We show that initially nonspherical 
seeds acquire spherical shape with passage of time. Applications to the description 
of the first-order phase transitions in nuclear systems, such as the nuclear gas-liquid 
transition occurring in low energy heavy-ion collisions and the hadron-quark tran- 
sition in the high energy heavy-ion collisions are discussed. In both cases we point 
out the important role played by effects of viscosity and surface tension. It is shown 
that fluctuations dissolve and grow as if the fluid were effectively very viscous. Even 
in the spinodal region seeds may grow slowly due to viscosity and critical slow- 
ing down. This prevents the enhancement of fluctuations in the near-critical region, 
which is frequently considered as a signal of the critical point in heavy-ion collisions. 



1 Introduction 



The description of first-order phase transitions is usually based on phenomeno- 
logical approaches. One constructs a thermodynamical potential depending on 
an order parameter, similar to that in the Landau theory of phase transitions. 
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Exploiting the fact that the time evolution of the collective mode is usually 
slow compared to the dynamics of microscopic processes one introduces a 
time-dependent equation for the order parameter pQ 

drf> = -r{A){6F/8</>), A it ■ .... ■ C (1) 

Here d is the dimensionality of the space, t is the time, F is a thermodynamic 
functional expressed in proper variables (e.g., the <fi- dependent correction to 
the Helmholtz free energy in T, V variables), T is the temperature, V is the 
volume and is an order parameter. For a non- conserving order parameter 
T(A) = cio and for a conserving one T(A) = a±A. Appropriate quantities 
should be chosen as order parameters for different physical systems under 
consideration. These could be a macroscopic wave function for metallic super- 
conductors and non-relativistic superffuids, as superfluid He 4 and He 3 , dipole 
moment for ferroelectrics, magnetic moment for ferromagnetics, a, u, p, it, 
i^-mean fields for various phenomena of nuclear physics, etc. 

There are many phenomena, where first-order phase transitions occur between 
phases with different densities. The description of such phenomena should be 
similar to that for a gas-liquid phase transition. Therefore it is worthwhile 
to find the solutions of corresponding hydrodynamical equations. Although 
some simplified analytical [2T3TH] and fragmentary two-dimensional numerical 
[5j solutions have been found, many problems still remain unsolved. In the 
general case one should construct an appropriate numerical scheme, which 
could describe the phenomenon. The hydrodynamical approach is fairly effi- 
cient for the description of heavy-ion collisions in a broad range of collision 
energy from SIS to RHIC energies (see e.g. [Dlr7f51l5irTU|ll|12lll3lll4fl5|16j ). As 
a simplification, most of three-dimensional hydrodynamical models use ideal 
hydrodynamics. Effects of viscosity and thermal conductivity are simulated 
then with the help of an artificially introduced friction between different com- 
ponents of the liquid. See, e.g. Ref. [7J for a three-fluid hydrodynamical scheme 
applicable in a broad energy range from SIS to SPS. 

There are arguments [21I3H] that the dynamics of a first-order phase transition 
is controlled by non- zero values of the kinetic coefficients. For nuclear systems 
transport coefficients are poorly known. Transport coefficients in nuclear mat- 
ter have been evaluated in [TTirT8] for temperatures and nucleon densities rel- 
evant for the nuclear gas-liquid (NGL) phase transition. Estimates [17] have 
demonstrated that the bulk viscosity is much smaller than the shear viscos- 
ity. Equations of non-relativistic non-ideal hydrodynamics have been solved to 
construct a description of heavy-ion collisions at SIS energies |6]. Some models, 
[TT][T2"] . describing the expansion of matter at RHIC energies solve equations 
of relativistic non-ideal hydrodynamics in two spatial dimensions and indicate 
that effects of viscosity in the state of strongly coupled quark-gluon plasma 
(sQGP) are minor. One concluded that the ratio of the shear viscosity to the 
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entropy density is r]/s < 0.2, see [12] . For certain materials, such as helium, 
nitrogen and water the experimental values for this ratio 77/s show a mini- 
mum at the phase transition [TH]. The properties of the bulk viscosity near 
the hadronic-to-sQGP phase transition critical point (endpoint) are discussed 
in Refs. [20|21|I22] . For different models have been derived results which dis- 
agree with each other. For example, Refs. [20] argue for an essential increase of 
the bulk viscosity at the critical point, contrary to the results of [21J. Lattice 
QCD calculations [23] applicable to the case of the sQGP demonstrate that 
the bulk viscosity is much smaller than the shear viscosity. 

It is expected that at large baryon densities and not too high temperatures the 
hadron-sQGP phase transition is of first-order pMf25] . Lattice results [261127] 
support this conclusion. Signatures of such a transition might manifest them- 
selves in heavy-ion collisions in a broad interval of energies, like those available 
at SPS (CERN) and at future low energy campaign RHIC (Brookhaven), and 
FAIR (GSI) and NICA (JINR) facilities. One expects strangeness trapping 
and enhancement of the kaon multiplicity fluctuations [28], enhancement of 
the soft pion yield [29] and baryon number density fluctuations [30], as signals 
of the first-order phase transition. The dynamics of classical modes of the chi- 
ral condensate in presence of an expanding fluid of quarks was studied in [TB] 
in terms of the cr-model. Quarks were treated in the framework of relativistic 
ideal hydrodynamics. It was indicated that significant density inhomogeneities 
appear around the critical point of the first-order phase transition. Ref. [25] 
argued that the fireball may linger longer in the vicinity of the critical point 
due to a divergence of susceptibilities, e.g., the specific heat. Refs. [2"9"lfT5] paid 
attention to the critical slowing down that limits the growth of the a-field 
correlation length in the vicinity of the critical point. Some models speculate 
about explosive freeze-out assuming an increase of the viscosity close to the 
critical point, see [3T] . 

Another relevant phenomenon is the NGL first-order phase transition that 
manifests itself in heavy-ion collisions in the expansion stage at low densities 
[32] . Possible effects of super-cooled vapor and super- heated liquid phases, as 
well as those of the spinodal region, have been considered in [33 34J. The occur- 
rence of a negative specific heat (Cp) was reported, as the first experimental 
evidence of the liquid-gas phase transition in heavy- ion collision reactions |35j . 
The complete list of references is too long to be included here. For a review 
of this interesting topic, see [36] . 

Mixed (so-called pasta) phases may occur in systems with two or more con- 
served charges (e.g. electric and baryon charges) undergoing first-order phase 
transitions, see [37]. The physical origin of this phenomenon is based on the 
fact that the electric charge can be conserved globally rather than locally. 
The pasta (gas-liquid) phases appear in the inner crusts of neutron star. Be- 
sides, pasta phase may arise also in the interior regions of compact stars, 
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provided the equation of state (EoS) allows for first-order phase transitions to 
the pion or the kaon condensate, or to quark matter [38]. In equilibrium, pasta 
is constructed of Wigner-Seitz cells. The conservation of the electric charge 
is assumed within each Wigner-Seitz cell. Each Wigner-Seitz cell contains a 
seed (nucleus) of one phase surrounded by matter of the another phase. Seeds 
are droplets for a small concentration of the new phase, and rods and slabs 
for higher concentrations. Only little is known how pasta phases are cooked 
dynamically, see [3T39] . 

In this paper, we will describe the dynamics of first-order phase transitions 
by means of the standard system of equations of non-ideal non-relativistic 
hydrodynamics: the Navier-Stokes equation, 



pd t Ui + p(uV)ui 

= -VjP + V fc ^j]{VkUi + ViU k - ~5j fe divu) + (5 ik divu^ , (2) 
the continuity equation, 

d t p + div(pw) = 0, (3) 
and the general equation for the heat transport, 



ds 

T dt + div ( s ") 

2 

= div(KVT) + r](y kUl + V iUk - -5 ik divu) 2 + C(divw) 2 . (4) 

a 

Here p = mn, m is the mass of the constituents, n is the density of the 
conserving charge (e.g, the baryon charge); P is the pressure; r] and ( are the 
first (shear) and second (bulk) viscosities; u is the velocity of the element of 
the fluid; S is the entropy, s = dS/dV; k is the thermal conductivity; as before, 
d is the dimensionality of space. We solve these equations numerically in two 
spatial dimensions, d = 2, and analytically for arbitrary d in the vicinity of 
the critical point. Then we perform estimations for the cases of the NGL and 
the hadron-sQGP phase transitions. 

The paper is organized as follows. In sect. [2] we analytically treat the dynam- 
ics of a system in the vicinity of the critical point of the first-order phase 
transition. First we perform a reduction of the general system of equations 
of non-ideal hydrodynamics to equations for the order parameters. Then, in 
order to solve the problem analytically, we use a density expansion of the Lan- 
dau free energy functional in the vicinity of the critical point and derive the 
equation of motion for the density variable in dimensionless units. We study 
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evolution of density fluctuations in metastable regions: the super-heated liquid 
and the super-cooled vapor. Then we consider the evolution of the seed shape 
and study the dynamics of density fluctuations in the spinodal region. In sect. 
E] we find numerical solutions of the general system of hydrodynamical equa- 
tions ([2]) - (j3J) in two spatial dimensions and compare these numerical results 
with those of the analytical treatment. We also pay attention to the specific 
patterns, which are not seen in a simplified analytical formulation, such as the 
dependence of the time evolution of the seed on its shape and on the power 
of the inflow of the surrounding matter. To be specific, throughout numerical 
calculations and analytical estimations we use phenomenological modified Van 
der Waals (mVW) EoS. Its expansion in the vicinity of the critical point is 
done in the Appendix A and implies validity of the mean field approximation. 
A modification of the VW EoS which we do does not change the universality 
class, being the same for many substances named the VW fluids. In Appendix 
B assuming thermal equilibrium we evaluate a fluctuation region. Although 
our consideration is very general, allowing for further applications of the re- 
sults to the description of specific nuclear dynamics, in sect. H]we specify the 
parameters first for a nuclear matter system undergoing a NGL phase tran- 
sition and then for a system undergoing a hadron-to-sQGP phase transition. 
Sect. formulates conclusions. 

Throughout the paper we use units h = c = 1. The results of this paper are 
briefly summarized in the letter |40j . 



2 Small overcriticality 

2. 1 Reduction of equations of nonideal hydrodynamics to equations for the 
order parameter 

Assume that EoS allows for a first-order phase transition and conditions are 
such that the system is somewhere in the vicinity of the critical point. In this 
case pressure isotherm as function of the density has a convex-concave shape, 
as for the Van der Waals EoS. 

Further in order to construct hydrodynamical description we need an expres- 
sion for a thermodynamical potential depending on appropriate thermody- 
namical variables in each space-time point. Working in (T, p) variables one 
may use relation between the pressure P and the Helmholtz free energy F: 

P = p(5[F(T,p)}/5p)\ T . (5) 
If one wanted to work in (s, p) variables, one could use that 
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P = p(5[E(8,p))/5p)\ s , T = p-\6[E(s,p)]/6s)\ p , (6) 

where E(S, V) is the energy. These expressions generalize standard thermo- 
dynamical relations P = -(dF(T,V)/dV)\ T = -(dE(S,V)/dV)\ s , T = 
(dE(S, V)/dS)\v to the slightly spatially inhomogeneous configurations. 

To treat the problem analytically let us expand quantities, entering EoS and 
equations of hydrodynamics, near some reference point (T r , p r ) on the curve 
P(T,p) (or (s r , p r ) on the curve P(s,p), respectively), where (T r ,p r ) are as- 
sumed to be close to the values in the critical point, i.e. < T cr — T r <C T cr , 
< (Puq ~ Pgas)/Pr "C 1. It is convenient to take p r satisfying the condi- 
tion (d 2 P(p,T)/dp 2 )\ T = 0, or, as an alternative, p r = |(p[f g c + pf^) can 
be chosen, where pj^ c and p^ are densities at the liquid-gas equilibrium 
(puq = Pgas = Pmc, P is the chemical potential) determined by the Maxwell 
construction (equal squares on the plot P(1/p)\t, cut off by the horizontal 
line P = const). Further to be specific we will take T r = T cr ,p r = p cr . Then 
the above mentioned condition (d 2 P(T, p)/dp 2 )\ r = is fulfilled. In case of a 
mVW EoS, which we exploit in this paper, all necessary explicit expressions 
are presented in Appendix A. 

We will consider evolution of fluctuations. These are, e.g., seeds of one phase 
in another phase and fluctuations like waves. Seeds can be of different shapes. 
Simplest forms of seeds in d — 3 spatial dimensions are spherical droplets and 
bubbles, rods and slabs. Also seeds can have more peculiar shapes. 

Since the evolution of collective modes is slower than that for microscopic 
modes, we further consider small velocities u(t,r) of the growth/damping of 
the density and temperature fluctuations. Thus we linearize hydrodynamical 
equations in the velocity "u" and in the density Sp = p — p r and temperature 
5T = T — T T variables. Introducing auxiliary variable z = div u and applying 
operator "div" to both sides of the Navier-Stokes equation we obtain 



p T (dz/dt) = -A[5P-(dr) T + ( T )z\, d = 2(d-l)/d. (7) 
Continuity equation becomes 



(d5p/dt) = —p T z. 
Replacement of z into Eq. ([7j) produces [4] 



d 2 5p 



A 



6[F{T > sp)] \ T+ c +P7 Hdv I+ ty sp 



5{5p) 



dt 



(9) 



where we expressed the pressure in terms of the free energy in T, p variables. 
We added an additional constant term C in square brackets which will be 
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specified below. Obviously AC = 0. Note that for Sp > (p > p cr ) we deal 
with one phase and for 5p < (p < p cr ), with another phase. Transport 
coefficients can differ in those phases. To simplify the problem we assume a 
smooth density and temperature dependence of transport coefficients at the 
transition through the critical point although in a narrow fluctuation region 
some of these quantities could have a singular behavior. We ignore this com- 
plication since such a temperature region is narrow and since it may take a 
long time to develop these singularities. 

Equation for the density should be supplemented with the equation for the 
entropy, which in case of slow evolution of the seed acquires the form 



T r (ds/dt) = c VjI (d5T/dt) = n r A5T, (10) 
where T] T , ( r and k t are the kinetic coefficients taken at T = T r and p = p r . 
In terms of s, p variables Eq. (jH]) reads 



d 2 5p 



A 



S[E(s,6p)] 
8{5p) 



+ p r 1 (dr] T + Cr 



85 p 
) ~dt 



and Eq. (TTDT) acquires the form of the equation for the conserving order pa- 
rameter, also known as Cahn-Hilliard equation [41] : 



T r (ds/dt) = K T A(5[E(s,p)]/5s)\ p . (12) 

We should note that Eq. differs from that is usually exploited in the 
framework of the phenomenological Landau approach, see (pQ), and from equa- 
tions used for the description of the dynamics of first-order phase transitions 
in heavy-ion collisions, e.g. see [I2 ; 43 : 44 : 29], and in relativistic astrophysical 
problems [IS]. Difference with Eq. ([T]) disappears, if one sets zero the square 
bracketed term in the r.h.s. of (jH]). Then Eq. becomes 



dSp S[F(T,6p)} S[F L (T,6p)] \ m ,n 



From the first glance, such a reduction procedure is legitimate, if space-time 
gradients are small. However for a seed, being prepared in a fluctuation at 
t — with a distribution 8p(t — 0, r) — 5p(0, r), the condition ®MM2 | t=Q ~ o 
should also be fulfilled. Otherwise there appears a positive kinetic energy con- 
tribution. Probability of such fluctuations should be suppressed. On the other 
hand, two initial conditions cannot be simultaneously fulfilled, if equation con- 
tains time derivatives of the first-order only. Thus there exists an initial stage 
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of the dynamics of phase transitions (t ^ tmit), which is not described by the 
standard Landau equation, see Eq. (TJP or < f73j) . being broadly exploited in con- 
densed matter physics. R.h.s. of (TT5]) presents thermodynamical force driving 
the system to the final (equilibrium) state. This force should become zero for 
t — > oo, i.e. when the system reaches the final (/) equilibrium state. Therefore 
C = —Pf- Thus instead of the Helmholtz free energy we may introduce the 
Landau free energy requiring that 



S[F L (T,5p)} 

Pr S(6p) \r- p - p f-° ( 14 ) 
in the final equilibrium state. 

Let the time scale for the relaxation of the density, following Eq. is t p and 
the time scale for the relaxation of the entropy/temperature, following ffTUj) , 
is tx- The latter quantity is estimated as 



t T = R 2 c v , r /K r . (15) 

Thus tx grows oc R 2 with increase of the size of the seed R. On the other hand, 
following Eq. (J9]) t p oc R (we show below that a seed of rather large size grows 
with constant velocity). Evolution of the seed is governed by the slowest mode. 
Thus, dynamics of seeds with sizes R < _Rf og (for tx(R) < t p (R)) is controlled 
by Eq. (Q for the density. Here i?f og is the typical seed size at which t p = t T - 
For seeds with sizes R > -Rf og , t T oc R 2 exceeds t p oc R and growth of seeds is 
slown down. Thereby, number of seeds with the size R ~ _Rf og grows with time. 
If conditions are such that k is sufficiently large, tx exceeds t p only for seeds of 
rather large sizes. At terrestrial conditions, when growing droplet in the cloud 
becomes sufficiently large and heavy, it falls down under the action of the 
gravity. Falling down, the seed causes an avalanche of secondary droplets, if 
there are accumulated already many droplets of the size R ~ -Rf og - Thus solving 
Eqs. |P)), n~0\) we are able to describe such a phenomenon as rain, provided 
gravitational forces are incorporated. Contrary, if k is sufficiently small, tx 
exceeds t p already for seeds of rather small sizes (at terrestrial conditions, it 
may occur when gravity is not yet efficient). Since for R > Rf og growth of 
seeds is slowing down and thus number of seeds with the size R ~ i?f og is 
increasing with time, there appears the fog. 

Note that seeds of the new phase are produced in the old phase owing to 
short-scale fluctuations. The latter fluctuations are not incorporated in above 
hydrodynamical equations describing by the mean field variables. Contribu- 
tions of short-scale fluctuations can be simulated by a random force induced 
in Eqs. ([9]), ({TO]) with the help of the 5-correlated source terms, cf. [2]. Being 
produced owing these source terms, large scale fluctuations (seeds) evolve in 
time following hydrodynamical equations. 
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2.2 EoS in the vicinity of the critical point and equation for the density in 
dimensionless units 



Consider evolution in <i-dimensional space of seeds in case when the heat 
transport is not yet efficient and the dynamics is controlled by Eq. for the 
density variable. For the sake of simplicity let us use a convenient parame- 
terization of the Landau free energy, i.e. the generating functional in Sp, ST 
variables, which variation in Sp produces equation of motion 



<)!■). = I — 

Pr 



\v 2 {S P f 



— eSp 



(16) 



where SFl = Fl[T, p] — Fl[T t , p T ], that allows for the first-order phase transi- 
tion. Here coefficients c > 0, A > are some functions of T r and p r , v and e 
are yet functions of 5T and |e| <C At> 3 . As we have mentioned, for convenience 
we choose T r = T cr , p r = p cr . For the mVW EoS, that we further exploit in 
our numerical calculations, see Appendix A, the value v 2 diminishes towards 
the critical point as v 2 oc (T cr — T) . 

For slightly inhomogeneous configurations the pressure can be expressed as 



SP = P-P(T I ,p T ) 

BP 1 B 2 P 

* Oflr^T + -^k Pr (5T) 2 - W5p + \{5pf - cA5p. (17) 



Here we used expansion of SP near T r , p r and Eqs. (fill) . (|T6|) . Thus 



e = P f - P(T r ,p r ) - ^\t^5T - -^k Pr (5T) 2 - ... 
-Pf - Pmc - n CT (p in - (if), 

Pi n and pf are the chemical potentials of the initial and final configurations (at 
fixed P and T). For the mVW EoS the pressure expansion in ST, Sn near T cr , 
n cr is performed in Appendix A. Then maximum value e max = P max — P uc ~ 
n cr (p max — Pmc) oc (T cr — T) 3 / 2 , where p max and the chemical potential p max 
correspond to the state where P(n) has local maximum p max ; and Pmc> A^mc 
are quantities on the Maxwell construction. 

The Landau free energy density ST re i = Fl(T CT) p cr ) and the value SP re i = 
P* 6(Sp) \t/P{T t , Pt), for T r = T cr , p T = p cr , for spatially homogeneous con- 
figurations constructed following Eq. (fl6|) . as functions of the density Sp, are 
schematically shown in Fig. [1] left and right, respectively. For e > (solid 
lines) the liquid state is stable and the gas state is metastable, and for e < 
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Fig. 1. The Landau free energy density 8T re x = 5Tl/ '^LiT cr , p cr ) and the value 
fiPrel = Pcr ^ F s[sp) P ^ It/ P(T cr , Per) for uniform matter, as functions of the order 
parameter 5p for the EoS determined by Eq. (|16|) . Dash horizontal line (e = 0) in 
the right panel shows Maxwell construction. 

(dash-dotted lines) the liquid state is metastable, whereas the gas state is sta- 
ble. For e = two minima of the Landau free energy coincide and correspond 
to the Maxwell construction on the curve 5P(l/p) (shown by horizontal lines 
in the plot 5P(5p) in the right panel). Using Eq. f|T7|) . we rewrite Eq. ([9]) as 



d 2 5p 

at 2 



cA5p + Xv 2 5p - X(5p) 3 + e - p r 1 (dr) T + Cr 



05 p 
~dt 



In dimensionless variables 5p = vip, & = Xi/l, i — 1, • • • , d, r — t/to, Eq. (JT 
is simplified as 



d 2 i) 

f 2c 

\\T 2 



A { (A^ + 2^(l-^ 2 ) + e 

2(*7r + Cr) ~_ 
AV z p r 



(19) 



1/2 



2e 



_ CPr 
(dVr + (r) 2 ' 



Eq. (jl~9l) is the key equation for our subsequent analysis. For p T = p cr , T T = T cr 
taken, as the reference point, and for the VW fluids (see Appendix A) pa- 
rameters demonstrate the following temperature dependence near the critical 
point: v oc \5T\ l l 2 , e oc (5Tfl 2 (e max = VSn cr \5T\ 3 / 2 /T^J 2 , e™ ax = 4/(3^)), 
and / oc |5T| -1 / 2 and to oc |5T| _1 . The latter value shows that relaxation pro- 
cesses in the vicinity of the phase transition critical point prove to be very 
slow. This means that close to the critical point it is, indeed, legitimate to de- 
scribe the phase transition dynamics in the framework of the non-relativistic 
approximation, as we do. 

Note that values c, A, v, e in Eq. (TT6T) are purely phenomenological coefficients. 
Introducing new variables A =Xm 2 ,v —v/m,r] =r)/m,( —(/m,e = e/m 
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one could exclude, e.g., dependence on the mass of the constituent m from 
all dynamical characteristics of the system. The value m is not well defined 
quantity in case of the quark system. In new variables pressure and the Landau 
free energy are proportional to m, whereas all parameters entering Eq. (1191) 
do not already depend on m. 

As we have mentioned, Ref. [22] expressed an opinion that, if at some inci- 
dent energy the trajectory passes the vicinity of the critical point, the system 
may linger longer in this region due to strong thermodynamical fluctuations 
resulting in the divergence of susceptibilities, that may reflect on observables. 
Contrary, we argue that fluctuation effects of the vicinity of the critical point 
in heavy-ion collisions can hardly be pronounced, since all relevant processes 
are proved to be frozen for 5T — ► 0, while the system passes this region during 
a finite time. 

As one can see from Eq. (fTPl . the dynamics of the phase transition is governed 
by single combination (drj + Q of shear and bulk viscosities, cf. [13]. This ob- 
servation might be interesting in connection with discussions of the particular 
effects of shear and bulk viscosities in sQGP, see [T2j . 

The Landau free energy (fl6j) . being expressed in dimensionless variables, be- 
comes 



SFr 



2 i f„l,2 i\2 



There exist homogeneous stationary solutions of Eq. ffT9]) : 



^~±l + e/4, |e|/4<l, (21) 
corresponding to the Landau free energy density 



tfi « (22) 

Thus for e > the upper sign solution describes stable liquid and the lower 
sign solution circumscribes super-cooled vapor, see Fig. [TJ For e < the lower 
sign solution describes stable gas and the upper sign solution circumscribes 
super-heated liquid. 
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2.3 Dynamics of seeds in metastable and in stable surroundings, fort p ^> tr 



2. 3. 1 Seed density profile 

In d = 3 space one deals with configurations of different symmetries, e.g., with 
spherically symmetric solutions (c? sol = 3, droplets/bubbles), axially symmet- 
ric solutions (d so i = 2, liquid and gas rods) and one-dimensional solutions 
(dsoi = 1, liquid and gas slabs). To describe mentioned configurations we 
search two-phase solution of Eq. f|T9|) in the form, see [2113114] , 

^=Ftanh[f-f (r)]+e/4, (23) 



£ = £i for kinks (walls separating metastable and stable phases), £ = J^f + £f, 



for rods, and £ = y£f + £f + £i> f° r droplets/bubbles. For e > upper sign 
solution describes evolution of droplets (or rods and kinks of liquid phase) in a 
metastable super-cooled vapor medium. The lower sign solution circumscribes 
evolution of bubbles (or rods and kinks of the gas phase) in a stable liquid 
medium. For e < the upper sign solution describes evolution of droplets in a 
stable gas medium, whereas the lower sign solution circumscribes then bubbles 
in a metastable super-heated liquid medium. From the kink solution for £o ^ 1 
(thin wall) one can easily construct the solution for the slab (— £o < £ < £o) : 

$ ~ Tsgn£ tanh[£ - sgn£ • £ (r)] + e/4. (24) 

The boundary layer has the length |£ — £o( r )| ~ 1- Outside this layer correc- 
tions to homogeneous solutions are exponentially small. Considering motion 
of the boundary for £o( r ) S> 1, we may put £ ~ £o( r ) i n P3"]) . (|2"41) . Then 
keeping only linear terms in e in Eq. ffTO]) we arrive at equation for £,o(t), 



±o e ~ ^ ZT- (25) 



2 dr 2 2 f (r) dr 



Volume and surface contributions to the Landau free energy 

Substituting fl24l) in fl20l) . subtracting infinite constant term and supposing, 
as we have used above £o( r ) S> 1, we obtain 



27r 3 / 2 A 3 - d - 1 Aw 4 / d - 1 



r(4oi/2)r(i + (3 - d sol )/2) Pl [ d sol ^ 3 



1 2 

+ 3 e so + o'So 



(26) 



2A is the diameter, height of cylinder and the length of the squared plate for 
d so i = 3, 2 and 1, respectively; T is the Euler T-function. The first term in (1261) 
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is the volume term and the second one is the surface contribution, 5FL )Sur f . At 
fixed volume in d = 3 space, the surface contribution for droplets/bubbles is 
smaller than for rods and slabs. Thereby, if a seed prepared in a fluctuation 
is initially non-spherical, it acquires spherical form with passage of time (see 
discussion in subsection 12.41 below). The surface term is 5Ft iSur f = aS, S is 
the surface of the seed, a is the surface tension. The gradient term in ffTBl) is 
as follows 

«VL = ^J(W - g/(Ve*)W« = = i«W. (27) 

Now we are able to express the surface tension through parameters of the EoS. 
For the mVW EoS, see Appendix A, we find 

a = o \5Tf\ a 2 = 32mn 2 cr T cr c- / = -— — , (28) 

bl cr n cr \ol \ L/Z 

5T=(T-T cr )/T cr . 

2. 3. 3 Role of the viscosity and surface tension 

There are only two dimensionless parameters in Eqs. (fT9l) . (I25I) . e and (3. 
Parameter e shows the difference in the Landau free energies of metastable and 
stable states, see Eqs. (122j) and ( 1231) . Dynamics is controlled by the parameter 
f3, which enters together with the second derivative in time. This parameter 
can be expressed in terms of the surface tension and the viscosity as 

(3 = (32T cr )- 1 [dr /r + Q- 2 alm. (29) 

The larger viscosity and the smaller surface tension, the effectively more vis- 
cous is the fluidity of seeds. For (3 1 one deals with the regime of effectively 
viscous fluid and at f3 ^> 1, with the regime of perfect fluid. Note that (3 does 
not explicitly depend on the n cr . As is shown experimentally [16], the ability 
of liquid domains to coalesce is controlled by an interplay between the surface 
tension and the viscosity. Parameter (3, which we introduced, is responsible 
for such an interplay. 

2.3.4 Critical radius 

For the case of a metastable seed prepared in a stable surrounding, both 
terms in the Landau free energy (|26l) are positive (at d so i 7^ 1). Thereby such 
seeds should shrink. For stable seeds developing in metastable surrounding the 
first (volume) term is negative. Therefore in this case there exists a critical 
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size of the seed, which can be found by minimization of (I26I) in £o : £(f = 
2(dsoi - l)/(3|e|). For the mVW EoS, see Appendix A, the condition ^ r > 1, 
that we exploited deriving Eq. (I25p . is fulfilled for |e| -C e 1 ™^. The closer initial 
state to the Maxwell construction, the smaller is |e|. 

The Landau free energy (1251) decreases with decrease of the droplet/bubble (or 
rod) size, provided £o(t) < £g r - Such seeds of the new phase, being produced, 
are dissolved with passage of time. For £o( r ) > £o r the Landau free energy 
decreases with increase of the seed size, that explains growth of overcritical 
droplets/bubbles and rods of the stable phase in the metastable medium. 

In dimensional units the critical size of the droplet /bubble or rod becomes 



i? cr (7e r 



^j^-yss (30) 

O"0 



2V3^n cr T cr \5T\ 1 / 2 ' 



where the latter equality is valid for the mVW EoS, see Appendix A, and we 
put |e| = '-ye max , where constant < 7 < 1 characterizes deviation of |e| from 
e max . The surface tension parameter o"o is given by Eq. fl28|) . The smaller |e|, 
the larger is the critical radius of the droplet /bubble (or the rod) growing into 
the new phase. For e — > (Maxwell construction) R cr — > 00. 

Slabs of the stable phase, being placed in a metastable medium, grow inde- 
pendently of what was the value of their initial size (R > I) . 



2.3.5 Probability of fluctuation 

The probability of the initial density fluctuation ip(r = 0) of the size £ ~ 
£o(t = 0) (in dimensionless units), see (124]) . is as follows 



W ~ e -^Ko(r=o)]/T ) (31) 

for £ ( r = 0) < £q t '. For £ ( r = 0) > £p r (unstable region) 5Fl decreases 
with the growth of however Eq. fl3Tl) does not hold anymore just indicating 
that these fluctuations may grow. Certainly, correctly calculated probability 
of appearance of a fluctuation with the size £ (t = 0) > ^ should decrease 
with increase of £o( r = 0). Being produced in random processes, seeds evolve 
then according Eq. (1251) . 

Two comments are in order. First, the profile (1241) describes only one density 
distribution among various possible configurations, which can be produced 
in fluctuations. We should consider fluctuations involving many particles with 
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the density p ~ pu q in the interior region and p ~ p gas in the exterior (provided 
liquid is stable phase). In difference with other configurations, the overcriti- 
cal droplet described by solution <HM gains in the Landau free energy and 
conserves the form of the boundary layer, tanh(x), with time. Fluctuations 
with a different shape of the boundary layer reach tanh(x)- like shape after 
passage of a time, ~ t rec , necessary for a reconstruction of the density pro- 
file. Second, various fluctuations may have distinct velocities ^ of the seed 
boundary, the latter being specified as the point, where ^£ = 0. However 
fluctuations containing a large number of particles, which we are interested 
in, can be cooked with not a negligible probability only provided the velocity 
^ is zero or very small (<^C |e|). Otherwise the Landau free energy in Eq. (|3T!) 
(being oc |e| -C 1, see ( !22l) ) would acquire essential positive contribution that 
would greatly diminish probability for the occurrence of such a configuration. 
Fluctuations with initially zero or very small values ^p, but with other density 
distributions in the surface layer compared to (T24"l) . will acquire a finite value 
of the velocity at the seed boundary after passage of a reconstruction time, 
~ tree; when the density profile reaches the form of the profile given by Eq. 
fflM . For t ^> t TCC we may simulate all these cases with the profile of Eq. ([Ml), 
assuming different values of ^ at r ~ r rec . 

2. 3. 6 Dynamics of slabs, d so \ = 1 

Consider dynamics of a slab of the stable phase (region < |£| < £o( r )) placed 
in a metastable surrounding (|£| > £o(t))- In this case one can find general 
solution of Eq. ([251) . First, let us obtain solutions satisfying initial conditions 
£o(t = 0) = £o(0), ^Ht=o — 0, since appearance of such seeds in fluctuations 
is more probable than, if ^| T =o were nonzero. Moreover ^| T =o = 0, if the 
seed is formed near the boundary of the system (provided the boundary is 
flat). Solution of Eq. ([25]) acquires the form 



Thus slabs of an arbitrary initial size (for £ 3> I) grow with passage of time. 
In dimensional units we obtain 



Ut) = U0) + 1\c\t-1 



e1/3[l-exp(-2r//3)]. 



(32) 



R(t) = R + u 



asymp 



3<V 2 \e\P 

2 1/2 A 3/2 V 4 



1 - 



exp 



( 



(drj T + C,)/3 



) 



(33) 



where 2Rq is an initial size of the slab and 




(34) 
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u asymp ( 1 e max )=l\ST\ 1 / 2 ^3/3T cr /(2m), 



Wasymp is the velocity of the growth of the slab at large values of time, t ^> £ init . 
Second equality (|34|) is valid for the mVW EoS. Then 



3lc l/2 e maxp ^ 

(35) 



\v 2 mn cr 72n cr T 2 r (dr] r + ( r )\5T\ 



(dVr + Cr)P a o m 

One can distinguish two stages of the evolution: an initial stage, t < t init 
ij ~ 7"init — P in dimensionless units), and a subsequent stage, t ^> i init (or 
t r init in dimensionless units). The initial stage lasts 



i«=%±^, (36) 
8/5(«i?7 r + Cr) (y\m 



9T cr n cr \5T\ Wn cr Tl(d Vr + Q\ST\ ' 
Second line (1361) is for the mVW EoS. For t <C tinit from ( 1331) we obtain 



i2(t) ~ R + wt 2 /2, w = 1 ' , (37) 

V2c 



that corresponds to the growing of the slab with constant acceleration. 

Another typical time scale is t p = -R/w a symp, since for t ^> t init from (1331) we 
obtain R = Rq + M a sym P ^ that corresponds to the growth of the size of the slab 
with constant velocity. In the vicinity of the critical point for the mVW EoS 
the velocity M a s y m P (7e ma:r ) oc \5T\ 1 / 2 and 



= r^^ R \ 5T \' 1/2 - (38) 

Thus time scales tj n i t oc and t p oc R\5T\~ l l 2 demonstrate that all 

processes freeze out at the critical point. 

For a large viscosity the time scale £ init oc l/(di] T + ( r ) is rather short (excluding 
the vicinity of the critical point) and the system rapidly reaches the asymptotic 
regime. Contrary, in case of an ideal liquid typical time scale t init is long and 
thus during a long time the size of the slab grows with acceleration. 
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For the initial condition corresponding to a nonzero velocity of the seed bound- 
ary (at r = T rec ) we find 



£o(t) = £o(r rcc ) + ^|e|(r - r rec ) - ^ ^|e| - ^| T=Trcc ^J 

x [l-exp(-2(r-r rec )//3)] (39) 

instead of (132]) . 

Another comment is in order. From Eq. (jSJ) using (|24|) and the condition that 
the hydrodynamical velocity (the velocity of the inflow of the surrounding 
matter) w in fl ow = deeply inside the slab, we find 



^inflowW = -ltd 1 (^ /rfr)[sg<tanh(^ - sgn£ • f ) + 1]. (40) 

At the right boundary of the slab u in Q OW (R = Ro(t)) = — It^^d^o/ dr) = 
-dR(t)/dt. For R > R (t) at large t, u inf i ow (R > Ro(t)) = -3eZt _1 = 
— 2it aflyreLp . In case, if the initial shape of the seed differs from that given by 
Eq. ([M]), the solution PHI) is invalid for t < t rec , where as before, t rec is the 
time scale of the reconstruction of the shape of the seed. From our numerical 
solutions below we will see that for effectively large viscosity, (3 -C 1, the 
reconstruction period proves to be rather short and in opposite limit, (5 1, 
of an effectively small viscosity this time interval is very long. 



2.3.7 Initial stage of the seed evolution (for arbitrary d so \) 

Let us return to the consideration of the general case of an arbitrary value d so i. 
As in previous subsection consider dynamics of stable seeds in a metastable 
surrounding. Then we may replace ±e to |e| in Eq. fT2o^) . We search solution 
of Eq. fTSoT) in the form of the Taylor expansion 

£o(r) = a + aiT + a 2 T 2 + a 3 T 3 + ..., (41) 
imposing initial conditions £o( r = 0) = £o(0), d ^| T - > | T= o = 0. Then we find 



£o(t) = 6(0) + l -w T r 2 (l - pj + .., r < min |r init , ^o(0)/w T } , (42) 



where w T = 2/3 1 



6>(o) 



For slabs this result follows from the general 



solution ( 1321) at t r init . In dimensional units Eq. ( 142]) is rewritten as 
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R(t) = R + ^(^l- JL^j + t < min {t init , (R /w) 1/2 } . (43) 
The acceleration 



w = d 2 R/dt 2 = (dsoi - l)Aw 2 (R - R cr ) (R cr R ) (44) 

changes sign at Rq = R cr . As we expected, seeds of undercritical size shrink 
with passage of time, whereas overcritical seeds grow. For seeds of a nearcritical 
size the process proceeds slowly (w oc \5T\(R — R cr )/R 2 cr ). For undercritical 
seeds of a small size, w oc —\5T\/R . Note also that w does not depend on the 
viscosity. For d so \ — > 1 from (143|) . (jUJ) we recover result (J3j 



One could think that for slabs in case of perfect fluid Eq. (14*3]) is valid for all 
times. Indeed, it follows from general solution (1331 for (3 — > oo (t; n i t — > oo in 
this case). But motion with constant acceleration, see (14*4"]) . becomes relativistic 
for large times. Deriving Eq. (IT"*"]) from general system of hydrodynamical 
equations we dropped quadratic terms in the velocity. Thereby deriving 



we additionally assumed that t <C yRo/w. 

Notice, if we supposed that the velocity of the seed boundary were finite at 
t = 0, we would obtain 



2 |C| ^o(O) rfr |r= ° 



dr ""7 (3 

instead of Eq. ( 1421) . compare Eqs. ( 145]) and ( 1391) for r rec = 0, r r ini t. 



r 2 + ... (45) 



^.5.5 Late stage of the evolution (t~^> t init ) of large seeds (R(t) ^> R cr ) 

Let (3 < Co - Then we may drop the term ^-4^- in Eq. (1251) and Eq. (I25p simplifies 

as 



dUr)/dr = - (d sol - l)^). (46) 

Note that solutions of Eq. (14(5]) do not satisfy necessary condition ^|r=o = 0, 
except for the case £o(0) = (,q t - Therefore initial stage r < r init should be in 
any case described by more general Eq. (125]) . which is of the second order in 
time derivatives (see also discussion above in subsection 12. II) . Thus besides the 
condition (3 -C £ 2 we should still require that r 3> r init = (3. 

For £o( r ) 3> £cT ( e -§- f° r initially overcritical droplets/bubbles and rods of a 
very large size, £o(0) Q r ) surface effects become unimportant. In this case 
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one can neglect the term oc l/£o in Eq. (146]) . Then for £o( r ) ^ £cf we arrive 
at the solution (1521) . Thus large seeds grow with constant velocity, 

3 

f (r) ~ f (0) + -|e|r, fl(t) ~ R + M asymp t, t > W (47) 
The time scale for the growth of the seed of the size R is 



r~r p (&) = ^, f ~ t p (fl) = (48) 

>J p | ^asymp 

for t p ^> tinit. Values t m it and t p are the same as in (136]) and (138]) . For the 
mVW EoS we estimate 



= ; 7^ = wfV^ri #- K T'JW 1 - (49) 

"asymp J vi cr -t cr 1 i | J~L cr 



Both t init and t p (R cr ^e max ) are oc I5T]- 1 . 



From ( 1361) and (1491) we see that for the system near the critical point the time 
scale t p (R cr ,e max )) is larger than £ init for (3 <C 1. In case of effectively small 
viscosity, (3 3> 1, the time scale t init exceeds the value t p (R cr , ■ye max ). Then the 
regime t p (R) 3> t init is reached only for R 3> R cr {^e max ). Note that inequality 
P < Ci is obviously fulfilled for all £ > £f for t p {R CT ^e max ) > £ init , i.e. 
r p(£cf) 3> /3, since the latter inequality can be rewritten as (3 <C (£o r ) 2 - 



,2.5.5 Evolution of seeds with nearcritical sizes (at t t- m i t ) 

Let c? so i 7^ 1 and we continue to consider evolution of seeds of a stable phase in 
a metastable surrounding. Consider a seed of initially nearcritical size, |£o(0) — 
£o1 ^ £cT\ a t an intermediate stage of its evolution, when the size of the seed 
still remains to be close to the critical one, |£o( r ) — £o(0)| "C Co"- We assume 
that conditions r 3> r ini t = (3 and (3 <C £q are again fulfilled. As we will see, 
evolution of the nearcritical seed proves to be very slow (cf. Eq. (1521) below) 
and there exists a time interval, where these conditions are fulfilled. Solution 
of ff46]) is then as follows 

t („\ i fcri £o( r ) — £cT /• M\ | (^sol — 1)T ■. 

?o(r) + 4o ln F7nl 7cf = &( ) + & ' ( 50 ) 

^o(r)/rfr = (deo(r)/dr)U = (4d - 1)(£o(t) - ^^oW] -1 . 

Using condition that the size of the seed is close to the critical size we find 
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Ur) = 6(0) + (4oi - W£o(0)]-%(0) - £?>■ (51) 
In dimensional units the latter equation renders 



DfA D , + (^ol - l)cp r (i? - Rcr) 

R{t) = R + U n . CI t, W n . C r = ~ ~ 7~ — • (52) 

RoRcr{dr] T + Cr) 

As from Eq. (|44|) . we see that seeds (droplets/bubbles and rods) of overcrit- 
ical size grow, whereas seeds of undercritical size shrink. The velocity of the 
growth/shrinking of nearcritical seeds (for d so \ 7^ 1, 16(0) — £o r |/£cT ^ 1) 
proves to be very low (oc \Rq — R cr \)- Typical time, when initially nearcritical 
seed still continues to be nearcritical, is 



, , (Rcr) 3 (dVr + C) ^ , 

t ~ tn.cr = 73 77 [35 5— f > tinit, {06) 

[dsol — L)Cp v \rCo — ttcr\ 

that justifies omitting of the term ^tt, as we have done it. In the vicinity of 
the critical point for the mVW EoS^ t n . cr ( 7 e maa: ) oc 7~ 3 |i? - R C r\~ l \5T\- z/2 ■ 
Thereby dynamics of seeds at this stage is very slow: t n . cr 3> (t p (R cr , r ye max ), t in i t ) 



2.3.10 Evolution of seeds of initially small size (I <C Ro <C R cr ) at t~^> t 



nut 



For r ri n it, ^0 ^ describing seeds of a small size (1 <C £o <C £o r , ^soi 7^ 1) 
we can drop the term oc e in (|46[) . Then solution satisfying initial condition 
£o(t = 0) = 6(0) acquires the form 



d£o(r) (4oi - 1) 



<*r v /6 2 (0)-2(rf sol -l)r 

In dimensional units 



R(t) ~ ^ - 2(4oi - l)tf 2 Ao- (55) 
For r > r init , £q(0) > /3, i.e. for t init <C i <C £dis> from fl55l) we find 



R(t) = R + Mdist, Mdis = -—77 T7- I 56 ) 

«o(ar? r + Cr) 

The time scale of the dissolution of the initial fluctuation of a small size 



(i _ (drjt + C)Rl _ l6n cr T cr (dr) T + ( r )Rl 
dls 2(d so i - l)cp r (4 i - l)al 
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(in dimensional units) proves to be oc Rq. Thereby, fluctuations of sufficiently 
small sizes are easily produced and then they are rapidly dissolved. As we see 
from (15"5"j) . (J2SD, at t > t init the speed of the dissolution first reaches a constant 
value and then the process is rapidly accelerated. 

Note that solving the problem we assumed £ 3> 1 and £q(0) ^> (3. Both 
conditions are satisfied at least in case of an effectively large viscosity, /3 -C 1. 



2.3.11 Numerical solution of equation for the droplet/bubble boundary 

Peculiarities of different regimes, which we have described analytically in limit 
cases, are demonstrated by numerical solutions of Eq. (1251) presented in Figs. 
|2] and [3] on example of the time dependence of dimensionless velocity ^ ■ 
Insertions in Figs. [2] and [3] show evolution of seeds at small times. 

Fig. [2] shows numerical solution of Eq. ( 1251) for undercritical seeds. In the 
left panel we consider evolution of the seed of a small size (£o(0) = 0.3£ cr ) 
and in the right panel, of the seed of a near-critical size ( £o(0) = 0.999^)- 
We see that undercritical seeds are shrinking with time. For the seed of a 
sufficiently small size (left panel) shrinking process is rather fast. The larger 
viscosity (smaller (3), the faster is the process. One distinguishes two stages of 
the process: an initial stage (r < r init = (3), and a dissociation stage, r init <C 
t ~ Tdis, see fl57j) . The latter stage is subdivided by two stages: r <C Tdi S 
(characterizing by approximately constant velocity) and r ~ t^s (when the 
process is rapidly accelerated). The curve labeled as "approx" demonstrates 
analytical solution Eq. (J51j) . being valid for r ^> r init , that requires r dis r init . 
The latter condition is not fulfilled for (3 = 10 (effectively small viscosity), 
but it is fulfilled for (3 = 10 _1 (effectively large viscosity). Deviation of the 
dash-dotted curve from the solid one for (3 = 10 _1 is due to the fact that the 
ratio £o(0)/£o r = 0-3 is not yet much less than unity, and £o(0) ~ 4 is not 
yet much larger than unity, whereas conditions £o(0)/£ C r "C 1, £o(0) 3> 1 are 
required for validity of the analytical solution (I54p . We have checked this with 
various choices of parameters. Significant difference between the dash-dotted 
and the dash curve in case of a small effective viscosity (large 0) is due to 
the fact that r init > Tdi S in this case. Thus for an effectively small viscosity 
dissolution of the small size seed occurs at the time scale r ~ r init . 

For the nearcritical seed (see right panel of Fig. [2]) one can distinguish three 
stages of the process: initial stage (r < r init = (3); an intermediate stage, 
when the seed still continues to have nearcritical size, (£g r — £o( r ))/£cT ^ 1 
(it lasts very long, up to r < 500); and a short subsequent dissociation stage. 
Dependence on the viscosity (the value (3) does not manifest itself at the 
intermediate stage. The "n.cr." arrow in Figure insertion shows the value of 
the velocity which follows from our analytical solution (j5ip . 
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Fig. 2. Numerical solution of Eq. (|25p for the undercritical droplet of initial size 
£o(0) = 0.3£g r (left panel) and for the near-critical seed, £o(0) = 0.999£o r , (right 
panel) for effectively large viscosity (solid curves) and for effectively small viscosity 
(dash curves); e = 0.1. Dash-dotted line is analytical solution given by Eq. ([53]) . 
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Fig. 3. Numerical solution of Eq. (125j) for the overcritical droplet of initial size 
£ (0) = 1.001Co r (left) and for f (0) = l.l£g r (right) for two values of j3. The arrow 
3e/2 shows asymptotic value of the velocity (in dimensionless units). 



In Fig. [3] we demonstrate numerical solution of Eq. ( 1251) for overcritical seeds. 
We see that these seeds are growing with time. For initially nearcritical seed 
(left panel), there exist four stages of the process. Initial stage lasts a time 
T ~ r init = P, see Figure insertion. Next is a prolonged stage, when (£o( r ) — 
£cT)/£(T <C 1 (velocity is almost constant for r < 200). Then there is a transi- 
tion stage (for r < 1500 2000), when (£ (r) - ^ r )/^ r ~ h and finally the 
asymptotic regime, when the velocity reaches new constant value |e (shown 
in Figure by arrow). On a so long time scale the memory about initial stage 
of the evolution (the latter depends on (3: r < r init = (3) is lost: dash and solid 
curves coincide. Note however that in dimensional units asymptotic regime is 
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achieved much faster in case of a more viscous fluid. 

The right panel of Fig. [3] shows evolution of the overcritical seed for £o(0) = 
l.l^Q r . The process proceeds faster, than for £o(0) = l-001^Q r , and in three 
stages. From initial stage through a transition stage the seed reaches asymp- 
totic regime. Transition regime lasts very long. 

As we see, presented numerical solutions of Eq. ff25l) support our simplified 
analytical considerations performed above. 



2.4 Evolution of the shape of seeds 



To be specific consider solutions in d = 3 space. Majority of seeds, being 
produced in fluctuations, have near spherical form. Then their shape evolves 
with passage of time. To describe this evolution let us expand £o{0,4>,t) in 
spherical functions: 



£o = ££o&m)*i,m(0,0). (58) 

l,m 

It is convenient to use shorthand notations £ (^ = 0, m) = t;o(l = 0) and £ (^ 7^ 

0, m) = £oz- Further assume that the seed has almost spherical form, i.e. £oz <^ 

1. Using that = + | J| + ^f 21 = l ( l+ ^ 01 ; in the linear approximation 
in £ i from Eq. (I25p we obtain equation for the time evolution of the seed 
boundary (for £ z)- ^ looks cumbersome. Dropping terms oc 0(^ a , ^pr) we 
are able to simplify equation as follows 



id% t _ (2-1(1 + im ol ld&i t>1 (59) 



2 dr 2 ^2(/ = 0) r) p dr 
P = P 



l + i/(/+l)£ - 2 (/ = 0,r) 



For £q(1 = 0, r) ^> I 2 , that we will further assume (just to deal with simpler 
expressions) one has /3 ~ (3. The value (,o(l = 0, r) obeys Eq. (J25l) . Undamped 
mode with I = 1 describes the motion of the system as a whole. 



2. 4-1 Initial stage of the evolution 

For r <C min (r init , \J(,o(0)/w T ^ , see solution (142!) . Eq. (159|) is further simplified 

as 
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1 R d 2 ioi , d£oi , (1(1 + 1) ~ 2)fo 

2^ + ^ + e 2 (/ = 0,r = 0) = °- (60) 

Solution of this equation is as follows: £oz = Cie XlT + C2e X2T . Using initial 
condition ^| r=0 = one obtains d = -^£ O z(0), C 2 = ^£o*(0), 



Ai,2 = ~\ ± 



I 2(2(1 + 1) - 2) (61) 



/3 V 2 W = 0, r = 0) 
Thus we find two damping solutions. For 



l<f (Z = 0,r = 0)<^2/3[Z(Z + l)-2], 1>1, (62) 

besides damping there occur oscillations. Inequalities (1621 are fulfilled for /3 > 
^q(/ = 0, r = 0) (limit of effectively very small viscosity). 



2.4-2 Limit of an effectively small viscosity 

The time scale of oscillations is r ~ r ^ ld / osc _ | m i n (^X 1)2 ) |— 1 ^ y/j3£ (l = 
0,t = 0) < r ini t provided inequality (|62|) is fulfilled. Condition Tz ld /° &c <^ Tjjdt 
is satisfied only for very large (3 ^> |e| -2 in case £o(^ = 0, r = 0) ~ ^ r and for 
(3 3> 1 for seeds of a small size £o(^ = 0, r = 0) > 1. In dimensional units the 
time scale 

3^J[l{l + 1) - 2] V T -I 5T I/ 
does not depend on the value of the viscosity. 

Damping occurs at the time scale T/ ld / damp ^ ^ that corresponds to 



t „ tf dam * ~ t init (64) 

in dimensional units. Thus the time scale for the reconstruction of the initial 
density profile of the seed (if initially it deviates only slightly from the profile 
Eq. (123]) ). is t rec ~ tf^ da,mp ~ t init in this case. 



2.4-3 Limit of an effectively large viscosity 

For /3 <C 1 (effectively large viscosity) the time scale min (r^it , 
\/£o(^ = 0, r = 0)/w T ) is very short and for r ^> r ini t one can drop the term 
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d ^ 2 °' ) in l.h.s. of Eq. (|59p . Then the latter equation is simplified as, cf. [2]: 



(2 - 1(1 + dfoi _ n ^ 

^ = o) dF" - (65) 

Its solution is as follows: 



fiW = ^(0)exp ||dr'[2 - /(/ + 1)]^ 2 (/ = 0,r ) j . (66) 

All modes with / > 1 prove to be damped. Thus an initially deformed seed 
acquires spherical shape with time. Typical time scale is 



t v = $(1-0,t-0) t n ^ P^l3 />i. ( 67 ) 

' [1(1 + I) -2] 1 p(d<r k + Q[l(l + l)-2\' 

We find that ~ tdis for not too high I and for £ (^ = 0, r = 0) < ^Jf, see Eq. 
( 157j) . The larger Z, the more rapid is the process. For £ (^ = 0, r = 0) > £g r , 
= 0, r = 0) - ^ r ~ ^ r we obtain # ~ t p . 



5 Dynamics of density fluctuations in spinodal region 



Assume that the system is driven to the spinodal region (region between the 
minimum and the maximum of 5P(6~p)\t in Fig. [1] (right). In the spinodal 
region even fluctuations of an infinitesimal amplitude and size may grow. To 
demonstrate this assume that the density is such that the free energy 8Fl is 
close to its maximum (5Fl ~ 0). Then we may linearize Eq. (1191) dropping ijj 3 
term. Setting 



^ = -| + Re{^ T+ ^}, (68) 

where ipo is an arbitrary but small real constant, we find two solutions of 
linearized Eq. (fT9l) . 



7^(A;) = (2(3)~ 1 (-k 2 ± y/k*{l-4P) + 8Pk 2 ). (69) 

Growing modes correspond to the choice of "+"-sign and k 2 < 2. Most rapidly 
growing mode is described by the maximum value of 7^(fc) and k = k m : 
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kl = 2^rv ^ km) = 2vh~i (70) 

Thus the time scale for the growth of seeds in aerosol (spinodal region) is 



Uev = t /^(k m ). (71) 

Increase of the wave amplitude is stopped, when the matter is exhausted in the 
regions, where ipoe" 1 ^ cos(/c£) < 0. It occurs for ^ e 7 * T ~ 1. After that, seeds 
(regions ipoe 1 ^ cos{k^) > with typical size l/jf), for ipoe 1 * 1 ' ~ 1, begin to 
expand and flow together, until the stable phase is cooked. With logarithmic 
accuracy the time scale of the formation of the aerosol-like state is estimated as 
t aer , following Eq. ( J7T1) . The well known phenomenon illustrating exponential 
growth of fluctuations in the spinodal region is the boiling of the champaign 
occurring, when one opens the bottle. 



2. 5. 1 Limit of an effectively large viscosity 

For effectively large viscosity (/? <C 1) we get two solutions: 

^\k) ~ 2 [l - k 2 /2 - (3k 2 (2/k 2 - l) 2 ] , ^\k) ~ -k 2 /(3, (72) 

for k 2 3> 4/?. For k < a/2, ^\k) describes growing mode. The mode 7^(fc) 
is damped. The most rapidly increasing mode corresponds to the momentum 
k m = a/2/3 1 / 4 (see (jTDjl ). and (k m ) — 2. The time scale characterizing 
growth of this mode is 

tl er ~ t /2 = (d Vl + C r )[A^Vr] -1 = + Q[Vn cr T CT \5T\}-\ (73) 

Second equality in (1751) is valid for the mVW EoS. The larger the viscosity 
and the smaller — 5T, the slower is the time evolution. Since k m = a/2/3 1 / 4 , 
the size scale of seeds is 

iC k l/(V2(3^). (74) 

Thus the size of seeds in aerosol increases with increase of the viscosity. For 
k 2 > 2 both modes are damped. 

2.5.2 Limit of an effectively small viscosity 

In case of effectively small viscosity {(3 3> 1) we get 
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7i ,(k) ~ ±kp/^\ - k*/2, (75) 
and the maximum value ( , yT ax ) corresponding to k m — 1 is determined by 



7 — (1) = lip. (76) 
As we see, the time scale characterizing the growing mode, 

4er ~ toft* = 2c 1/2 /(Xv 2 ) » t , (77) 

does not depend on the viscosity in this limit. For the mVW EoS one has 
v 2 oc \5T\ and t 1 ^ oc 1/|5T|. Evolution then is slow near the critical point. 
Conditions for the validity of the non-relativistic approximation, which we 
have used, are fulfilled. For t > t^ T , growth of modes becomes exponentially 
fast, see ( 1551) . Since k m — 1, the typical size of seeds in aerosol is 

-^aer — -^aer ^ ^aer' 

Modes with k 2 > 2 correspond to oscillations. Such fluctuations do not grow 
to the stable phase. 

2.6 Evolution of the density fluctuations in the system close to equilibrium 



Assuming that ip is close to its new equilibrium value, ip eq ~ ±1 + e/4, we put 
ip = ^eq + Sip in Eq. (Tl9l) and linearize the latter equation in Sip: 



-^^A { (A 5 ,,-4,,-^). (79) 
Setting 

5iP = Re{ip e^ T+l ^}, (80) 
where ipo is an arbitrary but small real constant, we find 

j^k) = (2/3)- 1 (-k 2 ± ^(l-4/3)-16/% 2 ) . (81) 

Dependence on the quantity e disappears from (1791) and ( 15T1) . 



27 



2. 6. 1 Limit of an effectively large viscosity 

In case of effectively large viscosity (/3 <C 1) and for k 2 > 16(5 there are only 
damped solutions. For k 2 3> 8/9: 



Fluctuations with large k rapidly dissolve with time. Existence of long living 
short-wave excitations is unlikely in the viscous medium. However there remain 
long-wave damped oscillations, for k 2 < 16/3. 

2.6.2 Limit of an effectively small viscosity 

In case of effectively small viscosity {(3 ^> 1) we get 



that corresponds to oscillating and slowly damped modes near the final equi- 
librium state, rather than to unstable modes. Since tovT? does not depend on 
the viscosity and t (3 — > oo for 77, ( — > 0, in case of the ideal fluid rapid oscilla- 
tions continue till the energy is transported to the surface of the system (the 
process is governed by the heat transport) or till the energy is radiated away 
in the course of direct reactions. Thus in case of effectively small viscosity the 
stable phase is covered by fine ripples during some rather long period of time. 

2. 7 The Reynolds number and turbulence 

The transition between laminar and turbulent flows occurs, when the Reynolds 
number 



exceeds the critical Reynolds number Re cr , which is very large 1000). Here 
u is the mean fluid velocity, A is the characteristic diameter of the body 
embedded in the fluid and p^ is the density of the fluid. In particular problem 
of the growth of the droplet /bubble of the stable phase developing in the 
metastable surrounding, u means the velocity of the growth of the seed, pfi is 
the density of the metastable phase (~ p cr in case of small overcriticality) and 
A = 2R. Supposing R = aR cr , a = const, for the typical radius of the seed, 
taking u from Eq. (1511) and replacing these values in (I54j) we find 



7^)^-^ + 4), ^\k) ~ -k 2 /(3. 



(82) 




(83) 



Re = uApfi/r] 



(84) 
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Re = —ot(3. 



(85) 



We see that for relevant values of a < 10 factor Re may reach critical value 
only for unrealistically large values of the parameter f3 (for a tiny effective 
viscosity). For values of (3 and 5T with which we are concerned, one has 
Re -C Re cr . 



3 Numerical integration of the system of hydrodynamical equa- 
tions in one and two spatial dimensions 

3.1 Setup 

To diminish computing time we consider solutions in d — 2 space. As is 
seen from analytical solutions presented above, qualitative description of the 
dynamics of the first-order phase transition remains similar for d = 3 and 
d — 2. E.g., Eq. (!26]) continues to hold for d — 2, provided one replaces 
SF^ [d BO i = 2;1]/(2L) to 5F^ d_2 ^ [d so i = 2;1]. Further, we assume validity 
of the isothermal approximation neglecting the heat transport effects. In this 
case evolution is governed by Eqs. (T5]) and ([3]). This means that in the present 
paper we will not simulate the late stage of the evolution of fluctuations, when 
a fog-like state is formed. 

To solve the problem we need to know EoS and transport coefficients. Since we 
will focus on demonstration of a qualitative behavior of the system undergoing 
the first-order phase transition, to avoid extra complications we will consider 
simplest case assuming that the viscosity does not depend on the density and 
the temperature. As we have mentioned in Introduction, for nuclear systems 
shear and bulk viscosities, as well as the heat conductivity, are poorly known 
and most probably the bulk viscosity is smaller than the shear one. Note that in 
analytical expressions, which we have derived, viscosities enter in combination 
(dr) r + Cr)- Thereby, and in order to diminish uncertainties we further put zero 
the bulk viscosity and vary the shear viscosity in broad limits. 

The Landau free energy and the pressure of the uniform matter behave as 
shown in Fig. HJ Note that in general case the Helmloltz free energy den- 
sity of the uniform matter does not produce two minima. Rather it fulfills 
the double-tangent construction (at values of densities corresponding to the 
Maxwell construction on the isotherm P(l/p)). E.g., it is so for the original 
VW EoS and for the EoS of the relativistic mean field Walecka model. We 
use the mVW EoS, P = f(T)Pvw, see Appendix A. Extra function f(T) 
is introduced for generality since temperature dependence of the purely VW 
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EoS is too simple to describe behavior of nuclear matter. At constant T for all 
f(T) hydrodynamical descriptions are self-similar. E.g., the description of the 
dynamics in the framework of the original VW EoS can be obtained with the 
help of the scaling t — > tJf(T) and (dr) r + ( r ) — > (dr) r + ( r )/J f(T), as it fol- 
lows from Eq. (Q. In our numerical calculations we specified f(T) in order the 
Helmholtz free energy had two minima, as well as the Landau free energy. Such 
a modification is quite not necessary since hydrodynamical equations enters 
only gradient of the pressure, which in both cases has the same form. Moreover 
we use f(T) — > 1 for T — ► T cr . In this case all the results valid in the vicinity 
of the critical point do not depend on / and are the same as for the original 
VW EoS. If we wanted to recover /-dependence for f(T — > T cr ) = f ^ 1 in 
the vicinity of the critical point, we could do it with help of the replacement 

k ->• V/o, P Pfo in E Q- (USD- 

In the present paper we consider dynamics of the phase transition at fixed 
pressure at the system boundary. In order to describe possible first-order phase 
transitions in heavy-ion collisions one still should take into account expansion 
of the fireball, see [37] . Moreover, one needs to use a more realistic EoS, e.g. 
that constructed in [3H]- These problems will be considered in forthcoming 
publications. 

The surface contribution to the pressure (the so called Laplace pressure) is 
taken into account in the low gradient approximation: 

p = p mVW - c V 2 p, (86) 

with p = mn. The coefficient "c" can be expressed in terms of the surface ten- 
sion, see (1281) . We examine both one- and two-space dimensional solutions. In 
general case no assumptions about cylindrical symmetry are used. To illustrate 
the dynamics of the overcritical and undercritical seeds we consider a fluctu- 
ation in infinite matter in both spatial directions x and y. Initial conditions 
correspond to a stable spot placed in the homogeneous metastable medium. 

3.2 Evolution of disks (d so \ = 2) in d = 2 spatial dimensions 

First consider dynamics of an initially static axial- symmetric seed with the 
conserved number density given by 

p(x, y; t = 0) = p out + {p in - p ou t)<d(Ro -r), r = \]x 2 + y 2 . (87) 

Such a seed would correspond to the rod in case d = 3. Densities p in and p out 
are taken to be those in homogeneous phases. We call these configurations 
liquid or gas disks in dependence, if pi n > p out or vise versa. 
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Fig. 4. Isotherm for the pressure as function of the density for mVW EoS, see Ap- 
pendix A, with initial and final configurations shown by dots (left column). Dash 
vertical line corresponds to the Maxwell construction (MC) on the curve P(l/p). 
In upper panel initial state relates to stable liquid phase disk in metastable su- 
per-cooled gas and on lower panel it corresponds to the stable gas phase disk in 
metastable super-heated liquid. Middle column demonstrates time evolution of the 
density profiles for the overcritical liquid disk (upper panel) and gas disk (lower 
panel). Numbers near curves (in L) are time snapshots; r = \/ x 2 + y 2 , \ST\ = 0.02, 
L = 30 fm. Right column, the same for initially undercritical liquid or gas disks. 

For d = 2, expansion of the mVW EoS near the critical point is valid only 
for \ST\ <C 1/7, see Appendix A, and analytical solution fl23l is applicable for 
\ST\ -C 1/16 (£o ^> 1)- Therefore we first take T very close to the critical tem- 
perature, in order one could quantitatively compare results of computing with 
analytical expressions. Results are presented in Fig. @] for T/T cr = 0.98. Pa- 
rameters of EoS are chosen, as for the NGL phase transition: T cr = 18.6 MeV, 
n cr /n sat = 0.42, n sat = 0.16 fm 3 is the nuclear saturation density. The configu- 
ration is computed for values of kinetic parameters chosen as r\ ~ 3.2 MeV/fm 2 
and (3 ~ 12.6 (effectively small viscosity). For a large viscosity general behav- 
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Fig. 5. The same as in Fig. [4] but for \6T\ = 0.15, L = 5 fm. 



ior of solutions remains the same but velocity of the seed evolution proves to be 
significantly smaller. Some peculiarities of the j3 dependence will be illustrated 
in next figures. In the upper panel of Fig. H]we demonstrate the time evolution 
of initially liquid disk and in the lower panel, of a gas disk. Left column demon- 
strates initial and final configurations on the curve P(p)- In the middle column 
we show dynamics of the initially overcritical seed Rq = 0.3 L > R cr ~ 0.2 L 
and in the right column, of the undercritical seed Rq = 0.1 L, in units of a 
relevant length scale L = 30 fm. The time snapshots are shown by numbers 
near curves in units L (i.e. 2 means At = 2L, etc.). We see from the middle 
column that in case R > R cr (in this example R ~ 1.5 R cr ) disks slowly grow 
with time. The initially selected distribution f|57|) acquires tanh-like shape (see 
(1231) ) for t > 20L ~ 600 fm. This is very large time (although in dimensionless 
units it corresponds to a rather short time scale r > 6, of the order of the 
duration of the initial stage, see the right panel of Fig. |3|. As we see from 
the right column of Fig. HI seeds of an initially small size Rq < R cr (in this 
example Rq ~ 0.5 R cr ) dissolve for t > 20 L ~ 600 fm. This value agrees 
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with the total time of the shrinking process described by the dash curve in 
the left panel of Fig. [2j The typical time scale characterizing the dynamics is 

t ~ tinit ^ ^dis- 

We have checked that the time evolution occurs in a line with above analytical 
consideration. The only difference is that in analytical treatment of the prob- 
lem with initial distribution of the tanh- form, we do not get values p < p ou t 
for droplets and p > p out for bubbles. As follows from our numerical solution, 
due to infall of the surrounding matter to the disk surface during the shape re- 
construction, the density decreases in the liquid disk neighborhood below the 
value of the density in the homogeneous metastable matter and it increases in 
the gas disk surrounding above the value of the density in the homogeneous 
metastable matter (see the middle column of Fig. H]). 

Note that it is unlikely to find the system in heavy-ion collisions in so narrow 
vicinity of the critical point as we considered, \5T\ = 0.02, since the typical 
time of the fireball expansion is much shorter than typical time of the evolution 
of fluctuations that we found. Therefore in Fig. [5]we also demonstrate the time 
evolution of disks for T/T cr = 0.85. For easier comparison with Fig. H] initial 
and final configurations are selected at approximately the same deviations 
(P — P m i n ) / (Pmax — Pnin)- In this Figure we take L = 5 fm, as the length unit. 
The system described by the VW EoS is already rather far from the critical 
point at T/T cr = 0.85. Indeed, for T/T CT > T comp /T cr = § ~ 0.844 there 
appears a region of negative pressures, that may cause some extra peculiarities 
in the processes under consideration. Different physical situations occurring for 
T < T comp were discussed in [33] . In this work we will avoid further discussion 
of the regime T < T comp . Although for T/T cr = 0.85 the dynamics looks 
qualitatively the same, as in case demonstrated by Fig. HI the critical radius 
proves to be significantly smaller (R cr ~ 1 fm instead of 15 fm in previous 
case) and the time scale characterizing the process is reduced. Therefore for 
T/Tcr = 0.85 seeds evolve much faster compared to the case T/T cr = 0.98. For 
overcritical discs the initially selected distribution (!H7|) acquires the tanh-like 
shape (see fl2]|) for t > (20 + 40)L = 100 + 200 fm. The typical time scale is 
t ~ ti n it, which in case \ST\ = 0.15 is much shorter than in case |<5T| = 0.02 
considered above. Initial disks of a small size (Rq ~ 0.5 R cr ) almost disappear 
for t > 10L = 50 fm. In this case t init ~ t^ is . All these values of time scales are 
larger or of the order of the time scale characterizing the fireball expansion in 
low energy heavy-ion collisions. 

We further varied parameters T cr , n cr , rj, (3 in broad limits (in the range 
relevant for the NGL phase transition) and checked that it does not change 
the qualitative picture presented in Figs. HJ [51 We demonstrate it in Fig. [6] 
showing the evolution of the disk surface with time for R Q > R cr . The disk 
surface is specified as the boundary, where density achieves the critical value 
(p = Per)- In the left panel results are presented for T/T cr = 0.98. The velocity 
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Fig. 6. Solution R(t) for the liquid disk boundary, p(R(t);t) = p cr , for several values 
of the viscosity. In left panel all parameters except viscosity are taken the same as 
in Fig. H (ST = 0.02, L = 30) and in right panel, as in Fig. (ST = 0.15, L = 5). 

of the growth of the seed (slop of the curve) increases with increase of (3. Large 
values of the dimensional time presented in Fig. [6j t ~ 50L, correspond to the 
value of the dimensionless time r ~ 5y//3. The latter value corresponds to the 
transition regime in the right panel of Fig. [3] (for all values of (3 presented in 
Fig. EJ). In this regime the velocity of the seed surface follows linear law. The 
asymptotic regime is achieved at much larger values of time. In the right panel 
of Fig. [6] results are presented for T/T cr = 0.85. For t ^> t rec the behavior is 
similar to that in the left panel. However at smaller values of time (t < t rec ) for 
large values of (3 (effectively small viscosity) there arise peculiarities. These 
peculiarities are associated with reconstruction of the initial density profile 
(1571) . occurring at t ^ t rcc . During this reconstruction period the typical radius 
of the seed may decrease. The dash curve demonstrates the same behavior, 
as is seen from the density profiles for time snapshots 2L, 10L and 20L in 
the middle column of the upper panel of Fig. If during the reconstruction 
process the size of the initially overcritical seed becomes smaller than the 
critical size, it causes subsequent shrinking of the seed. Namely this case is 
shown by the dash-dotted curve in the right panel of Fig. [61 To check this 
statement we increased initial size of the seed at given (3 and the behavior 
became similar to that shown by the dash curve. The larger (3 and \5T\, the 
higher is the seed velocity, cf. Eq. (|34|) . For fixed (3, evolution of the seed shown 
in the left panel (T/T cr = 0.98, L = 30 fm) is slower than that demonstrated 
in the right panel (T/T cr = 0.85, L = 5 fm). Thus peculiarities of the case 
of effectively low viscosity ((3 ^> 1), which are seen in the right panel of the 
figure, are due to a larger inertia than in cases presented in the left panel. 

In order to show how much the seed dynamics is sensitive to the choice of 
parameters of the EoS we take an another parameter choice T cr = 162 MeV, 
n/n sat = 1.3, relevant for the hadron-sQGP phase transition. Results are 
presented in Fig. [7J for configurations with approximately the same (P — 



34 




2 
1.5 



■e- i 

Cl 



0.5 
Q 



1 , 1 




1 1 1 1 1 


i ■ i ■ i 

y.25 
















^5 




■ / 




25/ I 
7 W / ■ 

/ / V Aioo _ 






— *f"""i i 




i.i.i 




i.i.i 



'0 0.1 0.2 0.5 1 1.5 0.5 1 1.5 
P/P r/L r/L 

cr 

Fig. 7. The same as in Fig.[5]but for T cr = 162 MeV, n/n sat = 1.3, r? ~ 45MeV/fm 2 
and for (3 = 0.2. 




1 0,2 0,4 0,6 0,8 

r/L 



Fig. 8. Evolution of bands (d so \ = 1) of initially large (left) and small (right) sizes, 
r = \x\. Parameters are taken the same as in Fig. [5j 
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Fig. 9. The same as in Fig. El but for \ST\ = 0.15, L = 5 fm, n ~ 23 MeV/fm 2 and 
(3 ~ 0.042. 

Pm in) / (-Pmax — Pm in ) as shown in previous figures. As in Fig. [5] we take T/T cr = 
0.85 and compute the configuration for rj ~ 45MeV/fm 2 and for (3 = 0.2 (ef- 
fectively large viscosity). As we see, typical time scales and the shapes of 
configurations look similar to the case presented in Fig. [5l in spite of the 
parameter sets are completely different. 

3.3 Evolution of bands (d so \ = 1) in d = 2 spatial dimensions 

In Fig. [8] we show dynamics of liquid bands in metastable gas phase. Values 
\ST\ = 0.02, L = 30 fm, are taken the same as in Fig. [5] and we again choose 
T cr = 18.6 MeV, n cr /n sat = 0.42, rj ~ 3.2 MeV/fm 2 and (3 ~ 12.6 (effectively 
small viscosity). These solutions are similar to slabs in d = 3. Left panel 
shows time evolution of a band of a large initial size (Rq = 0.3L), whereas 
right panel demonstrates evolution of a band having initially rather small size 
(Ro = 0.1L). In difference with disks (solutions with d so \ ^ 2) in both cases 
(for large and small initial sizes of bands) dynamics looks similar: bands of the 
stable phase, being prepared in the metastable phase, undergo growth to the 
new phase. Nevertheless, we also see that during the shape reconstruction the 
slab first begins to dissolve and then grows. This peculiarity appeared since 
initial form of the density distribution that we exploit in numerical calculations 
deviates from the form given by analytical solution (124"|) . Thus even for slabs 
there might exist a small critical size, that depends on peculiarities of the initial 
density profile. Slabs having sizes smaller than this critical size could then 
completely dissolve. 

Fig. shows the same as Fig. [81 but for \ST\ = 0.15, L = 5 fm, r] ~ 
23 MeV/fm 2 and (3 ~ 0.042 (effectively large viscosity). We see that the qual- 
itative picture of the time evolution remains the same. 

The probability to prepare a band in a fluctuation is tiny. However, bands 
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of the stable phase could be formed near the system boundary, provided the 
latter is flat. 



In Fig. [TU] we demonstrate the law for the growing with time of the band 
boundary R(t) (in the left panel) and the velocity of the boundary u = dR/dt 
(in the right panel) for different values of the viscosity. As for discs, the band 
boundary is specified as the point, where the density achieves the critical value 
(p = Per)- Results are presented for T/T cr = 0.98, L = 30 fm. For small values 
of time (see Figure insertion) R(t) obeys the quadratic law, as it follows from 
Eq. fl32|) . Solid curve (effectively large viscosity, (3 — 0.1) follows the law fl32l) 
for t > t rec ~ 100 fm. The higher (3, the larger is deviation from this law since 
the reconstruction time increases then as t rec oc It is clearly demonstrated 
in the right panel, where we present the time dependence of the velocity of 
the seed growth. As follows from the Figure, even for large times the velocity 
u does not obey the scaling law, u oc y/]3 (as it would follow from Eq. ([3*3"]) ). 
This is so, because values of time t ~ 200L still correspond to the transition 
regime in the right panel of Fig. [3] (for all values of f3 presented in Fig. [TO]) . 
In this regime the velocity of the seed surface still slowly increases with time. 
The asymptotic regime is achieved at larger values of time (or for smaller /3 
at values of time shown in Figure). 



Another important issue is presence of the damped long-wave oscillations 
which are clearly seen for all values of the effective viscosity. They occur at 
t < t Tec . Besides, in case of effectively small viscosity short-wave oscillations 
are clearly seen. As follows from Eq. (!83l) . stable phase is, indeed, covered by 
fine ripples. 
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Viscosity/ Wave number 


Small 


Large 




oscillation 


damping 




growth 


growth 



Table 1 

Different scenarios of the evolution of initial disturbances in spinodal region. 
3.4 Evolution of fluctuations in spinodal region 

As the initial density profile, we take now a wave developing on the background 
of a density p ( the latter value is chosen somewhere in the spinodal region), 
cf. sect. |2j The time dependent solution has the form 

p(t)=p + A f(t)sm(kr), (88) 
where A is a small constant and /(0) = 1. 

We checked that in case of an effectively low viscosity (/? > 1) for k > k cr = 
\/2 there appear oscillating modes and at an effectively high viscosity 1) 
there are only damped modes. For k < k CT = \pl there are growing modes. 
Various regimes, as corresponding to growing, oscillating and damping initial 
disturbances, are indicated in Table 1. 

In Fig. [II] we show time evolution of wave amplitudes given by Eq. (l88~l) for an 
undercritical value of the wave number k (left panel) and for an overcritical 
value (right panel). We take T cr = 18.6 MeV, n cr /n sat = 0.42, c ~ 5.56 • KT 3 
fm 2 . In case of the overcritical value k and an effectively small viscosity (J3 = 
10) we demonstrate the change of the amplitude in the half-period of the 
oscillation. Such a behavior fully agrees with that follows from our analytical 
treatment of the problem, see ( 1721 . ( 1521 . ( 1531) . We have checked that in case 
of a small overcriticality (for \ST\ = 0.02) slopes of the curves coincide up 
to the third digit with values of 7 calculated in subsect. 12.51 Even in case of 
sufficiently large deviation from the critical point (for \5T\ = 0.15, as presented 
in Fig. [TTj) . the difference of the curves f(t) obtained numerically from those 
calculated analytically is less than 30%. 

In Fig. [T2] we show time evolution of the wave amplitudes given by (1551) . for an 
undercritical value of the wave number k (left panel) and for an overcritical 
value (right panel) for the parameter choice T cr = 162 MeV, n/n sat = 1.3, as 
for the hadron-sQGP phase transition. 

We see that in case relevant for the hadron-quark phase transition the evolu- 
tion is more rapid compared to the example of the configuration presented in 
Fig. [Til relevant for the NGL phase transition. Nevertheless even in the former 
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Fig. 11. Time evolution of the wave amplitudes f(t), see Eq. (|88|) . normalized to 
the amplitude of the initial disturbance. Solid line is for effectively small viscosity 
(f3 = 10) and dash line, for the large viscosity ((3 = 0.1). Left panel: the undercritical 
wave number k = 21 /L (growing modes). Right panel: the overcritical value k = 81 /L 
(oscillation modes for large (3 and damped modes for small 0). Other parameters 
are taken the same, as in Fig0 




case the time evolution remains sufficiently slow, especially at an effectively 
large viscosity. Then the typical time scale t ~ 10 fm is comparable with the 
total time of the fireball expansion. Thus in heavy-ion collisions during ex- 
pansion of the fireball the system may linger in the old phase for a while even 
at T < T cr . This means that the equilibrium value of the critical temperature 
of the phase transition might be significantly higher than the value which may 
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Fig. 13. Isolines of the density n/n cr demonstrating the time evolution of initially 
asymmetric disk of overcritical size (in both directions) for effectively large viscosity 
{(3 = 0.1). Density isolines are drawn with an increment equal 0.25. Time snapshots 
are indicated in the sub-figures, \5T\ =0.15, L = 5 fm. 



manifest in the growth of fluctuations in experiments. 



3.5 Evolution of asymmetric spots (d so \ = 2) in d = 2 space 



Now let us demonstrate how seeds having shapes significantly different from 
discs become the discs with passage of time. In Figs. [13] and [TH we show 
time evolution of the initially rectangular seed for effectively large ((3 = 0.1) 
and very small (/3 = 10 3 ) values of the viscosity, respectively. In case of an 
effectively large viscosity the shape of the seed monotonously transforms to 
the spherical one for typical time t rcc ~ 30 fm. For an effectively small viscosity 
the dynamics is a more peculiar. In a line with our findings, see Eq. ( |62|) . one 
can recognize oscillations of the form. In the process of oscillations some pieces 
of matter first decouple with the growing seed and then fly away. Besides, the 
process lasts longer than in case of a large viscosity, t rec ~ 150 fm. 
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Fig. 14. The same as in Fig. \13\ but for effectively very small viscosity (3 = 10 3 . 
4 The NGL and hadron — sQGP first-order phase transitions 

4-1 The NGL phase transition 

For the description of the NGL first-order phase transition we take values 
T cr ~ 18.6 MeV, n cr /n sat ~ 0.42, n sat ~ 0.16 fm -3 , similar to those one 
uses in the RMF models describing this transition, e.g. see [33]. We use the 
mVW model for the EoS, see Appendix A. Parameters of the EoS are then 
as follows: a ~ 3.14 • 10 2 MeV • fm 3 , b ~ 5 fm 3 , A ~ 2.85 • 10~ 6 fm 6 /MeV 2 , 
v 2 ~ 1.56 • 10 4 |<5T|MeV 2 /fm 6 , e max ~ 2.16(<5T) 3 / 2 MeV/fm 3 . 

In order to estimate the gradient coefficient "c" we suppose that l(T = 0) = 
ddif — 0.5 fm, see Eq. ( IT91) . where d^f is the diffusion length, as it follows from 
the Woods-Saxon parameterization of the density profile of the nucleus. Then 
c ~ 5.56- 10~ 3 fm 2 , and for an appropriate value \ST\ ~ 0.15 we find / ~ 1.3 fm. 
In order to obtain l(T = 0) ~ 0.5 fm we should take o"o — 4 MeV/fm 2 . Then 
for T = 0.85 T cr we get a ~ | MeV/fm 2 . Using these values and taking 
r] ~ 23 MeV/fm 2 following Ref. [18J, we obtain /?ngl — 0.023 that corresponds 
to the limit of effectively very large viscosity. 

Then we evaluate the time scale to — 20|5T| _1 fm, see Eq. ( IT9|) . Also we 
are able to estimate values t p , tdis as they follow from Eqs. (|48j) . (!57|) . We find 

t p (-fe max ) ~40i?7- 1 |<5T|- 1 / 2 ,t p ( J R cr .,7e maa: ) ~ 307- 2 |5T|- 1 , t dis ~ 20R(R/im). 
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The initial stage of the process occurring for t £ t- m it — 0.06|<5T| 1 , see (136]) . 
proves to be very short (i init <C t p ). 

Further following [18] we evaluate thermal conductivity k ~ 0.08fm~ 2 and 
estimate typical time for the thermal transport tr- From (fT5l) we find tx — 
3R(R/fm). Here we used that for symmetric nuclear matter cy ~ 3n for 
T 3> e-p and cy ~ 2(7r/3) 2//3 mn 1 / 3 T for T <C ep. Thus the heat transport 
might become operative for R > Ri og ( r ye rna ) ~ 137~ 1 |5T|~ 1 / 2 fm, where Rf og 
is a typical radius of seeds in the "nuclear fog". Using that R cr {^e max ) ~ 
0.97 _1 |<5T| _1 / 2 fm, see Eq. f l30l . we obtain R{ og 3> -R cr - Thereby there is a 
long time interval, where the heat transport is not yet efficient and solutions 
presented in this paper are valid. 

Assuming that the fireball reaches spinodal region, following f!73|) we may 
estimate the time scale for the formation of the aerosol, t^ er ~ 10|5T| _1 fm, 
and the size scale for seeds in aerosol, R*l er ~ 2 _1 /3 _1//4 Z ~ 0.6|5T|~ 1//2 fm. 
With the help of Eq. (1851) for (3 ~ 0.023, we evaluate typical Reynolds numbers 
Re ~ (0.1 4- 1) < Re cr ), for a ~ 1 4 10. 

We see that all relevant time scales for the formation of seeds are long (> lOfm) 
increasing, when the system comes closer to the critical point. Thus it is likely 
that only an initial and might be an intermediate stage of the NGL first-order 
phase transition may manifest itself in the course of heavy-ion collisions. Also 
it seems unlikely to observe effects of a narrow vicinity of the critical point of 
the phase transition, since the fireball is located in this region a short time com- 
pared with the time scale characterizing the growth of seeds. This also means 
that the critical temperature of the first-order phase transition calculated at 
assumption of the thermal equilibrium might be significantly higher than the 
value, which may manifest in the growth of fluctuations in experiments. Note 
that values of the critical temperature calculated in different models are in the 
range T cr ~ 15 4 20 MeV, whereas experimental value, as it follows from the 
analysis of the multi- fragmentation is T cr ~ 5 4 8 MeV, see [32] ■ One usually 
associates this difference with the finite size effects. We point out that at least 
partially the difference can be explained by the dynamical effects. 



4-2 The hadron-sQGP phase transition 

In case of the hadron-sQGP first-order phase transition critical values T cr 
and rather unknown and can be varied in a broad range. Values 

T cr = (80 4- 170) MeV and n cr = (1 4 6)n sat are used in different models. 
For rough estimates we take values T cr ~ 162 MeV, n cr /n sa t — 1-3, as they 
are obtained following lattice calculations, cf. [27]. We again use the mVW 
model for the EoS, see Appendix A. Parameters of the EoS are then as fol- 
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lows: a ~ 8.76 ■ 10 2 MeV • fm 3 , b ~ 1.60 fm 3 , A ~ 7.80 ■ 10- 5 <r 3 fm 6 /MeV 2 , 
v 2 ~ 1.56 • 10VI^|MeV 2 /fm 6 , e max ~ 58.4 • |5T| 3 / 2 MeV/fm 3 , m g is the ef- 
fective quark mass, g = m g /(300MeV). Further, we obtain l(T = 0) ~ 0.2 fm 
(radius of confinement) for uo — 40 MeV/fm 2 . If one used o"o — 100 MeV/fm 2 , 
one would estimate l(T = 0) ~ 0.5 fm. 

Next we estimate s(T) ~ 7T 3 (T/T cr ), c v ~ 28T 3 (T/T cr ) at T near T cr , as 
it follows from the lattice data [49J. Assuming minimal value of the viscosity 
Vmm = s/(47r) = 60MeV/fm 2 , ( min = we evaluate maximum value of (3: 
/3^qq P ~ 0.015g for a ~ 40 MeV/fm 2 [5U], that corresponds to the limit of 
effectively very large viscosity. Even for o~o ~ 100 MeV/fm 2 , m q = 600 MeV 
we would get (3^qgp — 0-2 "C 1- Note that following [20] the bulk viscosity 
diverges in the critical point. If were so (/3 — > 0), the quark-hadron system 
would behave as absolutely viscous fluid, like glass, in near critical region. 
Contrary, Refs. [2~l~lfl~2"] argue for a smooth behavior of the bulk viscosity. 

With/3 = 0.015, we further estimate t ~ 2|<5T|- 1 fm, t p (7e max ) ~ g.l^ -1 ? 1 / 2 ^^- 1 / 2 
and tdis — 14g(-Ro/f m )-Ro- The time scale for the formation of the aerosol is 
^acr — l^l -1 fr 11 ' an d the size scale for seeds in aerosol is R2 er — 0.24I5TI -1 / 2 fm. 
Only ti nit ~ 0.03g|5T| _1 fm proves to be small (excluding quite small ST). 

For the thermal conductivity we use an estimation k ~ a^/m, see [T7j. Fac- 
tor «o depends on the EoS used. We will take «o = 3. Then one recovers 
appropriate relation between values of k and r\ for NGL transition, which we 
have used above, see [TB]. For the hadron-sQGP transition this estimation 
renders k s qgp — 3r)/m q . Then we are able to evaluate the scale of the heat 
transport time, tr — 26g(i?/fm) 2 fm. The heat transport becomes operative 
for R > R fog (je max ) ~ O^^q- 1 / 2 ^^ 2 fm. Here R iog is the scale of size 
of the seed in the quark (or hadron) fog-like state. Using that R cr ( , y€ rnax ) ~ 
0.37- 1 |5r|- 1 / 2 ((Jo/(40MeV/fm 2 )) fm, for a = (40 100) MeV/fm 2 we ob- 
tain Rf og < R cr . Thereby, the heat transport might be always operative for the 
description of the evolution of overcritical seeds in hadron-sQGP phase tran- 
sition. The value i?f og proved to be very small (< 0.1 4- 1 fm). However typical 
time tj< is rather long. Therefore, the system most probably would have no 
time to fully develop a fog-like state in a hadron-quark phase transition in 
heavy-ion collisions. 

For the system located in the vicinity of the critical point all estimated time 
scales (except t init ) are very large. If the system trajectory paths rather far 
from the critical point, all time scales, except tx, become of the order or 
less than the typical life-time of the fireball (~ 10 fm at RHIC conditions). 
Reynolds numbers are < 1, being much smaller than Re cr > 1000. Thereby, 
turbulence regime is not reached. 
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5 Conclusion 



In this paper we studied the dynamics of systems undergoing first-order phase 
transitions. We formulated analytical description of the problem (in general 
case in d + 1-dimensional space-time) valid for systems in the vicinity of the 
critical point. The analytical solutions were derived for configurations of cer- 
tain symmetries (droplets/bubbles, rods and slabs for d — 3). 

Then, the general system of equations of non-ideal non-relativistic hydrody- 
namics was numerically solved for a modified Van der Waals equation of state. 
Results for the original Van der Waals equation of state can be obtained by 
simple re-scaling of the time and the viscosity. Since there exist many differ- 
ent regimes for the relevant processes, in this paper we partially restricted our 
analysis. We did not incorporate the heat transport, which governs evolution 
of seeds of a large size. Also we focused on the case, when concentration of 
seeds of the new phase in the old phase is still rather small and one can ignore 
their coalescence. For simplicity we did not generate fluctuations in random 
processes considering evolution of the given seed after it has been produced 
in a fluctuation. Generalizations will be presented elsewhere. With mentioned 
reservations we revealed and studied general features of the dynamics of first- 
order phase transitions. These main features are as follows: 

(i) Essential role in the dynamics of the first-order phase transitions is played 
by the viscosity effects. This might be very important, since existing three- 
dimensional hydrodynamical schemes pretending to study phase transitions 
exploit equations of ideal hydrodynamics, whereas viscosity effects are sim- 
ulated only implicitly, e.g., with the help of phenomenological coefficients 
responsible for a friction of fluids. 

(ii) Because of the surface tension (except for one-dimensional slabs) there exists 
a critical size for the seeds of the stable phase. Seeds of overcritical sizes grow 
with time, while undercritical fluctuations dissolve. The closer to the critical 
point, the slower are processes. Even far from the critical point overcritical 
seeds grow slowly. 

(iii) We have shown that seeds with an asymmetric shape become spherical 
with time. This process is slow. For systems with effectively large viscosity 
the shape of the seed changes steadily to the spherical one. For systems 
with effectively small viscosity the seed undergoes long-wave and short-wave 
damped oscillations in the process of acquiring spherical quasi-equilibrium 
shape. Short-wave oscillations damp very slowly. 

(iv) In the spinodal region the system is unstable against generation of waves 
(with not too high wave-number). However instability develops slowly, if the 
system is close to the critical point (the time scale tends to infinity at the 
critical point). Far from the critical point processes become more rapid. The 
latter observation can be very important for studying possible signatures of 
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the gas-liquid and hadron-quark first-order phase transitions in heavy-ion 
collisions. 

As a signal of the phase transition, some models suggest to search anomalies in 
the behavior of the derivatives of thermodynamic quantities, e.g. specific heat. 
These anomalies appear due to fluctuations, whereas we have shown that those 
anomalies in fluctuations may not have sufficient time to develop. Moreover 
we conclude that in heavy-ion collisions the system may linger in the old phase 
( e.g. in the QGP state, or in the gas state) longer during the fireball expansion, 
even when T(t) has already decreased below the corresponding value of the 
equilibrium critical temperature of the phase transition. In another words, this 
means that the value of the critical temperature calculated within equilibrium 
thermodynamical models might be significantly higher than the value, which 
may manifest in the growth of fluctuations in experiments. 

Typical dimensionless parameter f3 that separates effectively viscous and per- 
fect fluid regimes proves to be oc o"o/(§?7 + C) 2 > where a is the surface tension 
at T = and rj and ( are the shear and bulk viscosities. According to our esti- 
mates the system undergoing nuclear gas-liquid phase transition in the course 
of heavy-ion collisions at low energies represents effectively very viscous fluid 
(/? <C 1). There exist arguments that strongly coupled quark-gluon plasma 
state, which is, as commonly expected, formed in heavy-ion conditions at 
RHIC, represents almost perfect fluid in the cross-over region (see e.g. 
Estimating the ratio of the viscosity to the entropy density as rj/s < 0.2 it was 
concluded [12] that strongly coupled quark-gluon plasma is the most perfect 
liquid known. It is usually believed that this property of the plasma will sur- 
vive at finite baryon density for systems in the vicinity of the critical point of 
the first-order phase transition (critical end point). In contrast our estimates 
show that the system undergoing the hadron-quark first- order phase transition 
in the course of violent heavy-ion collisions represents effectively very viscous 
fluid ((3 < 1). 

We found that the heat transport effects may play important role in description 
of the hadron-s QGP phase transition dynamics, whereas these effects are much 
less pronounced in the case of the nuclear gas-liquid transition. 

In future we plan to use a realistic equation of state to study the heavy-ion 
collision dynamics. 

After our paper has been submitted to the journal there appeared interesting 
paper [51] devoted to the description of fluctuations in the spinodal region 
at the hadron-sQGP first-order phase transition which well complements our 
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study. 
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6 Appendix A. Modified Van der Waals EoS 



The best known example to illustrate principal features of the first-order phase 
transition is the Van der Waals fluid. The pressure is given by 



NT N 2 a nT „ 
' ^ = y^Tpjb ~ ~yT = ~\~^bn ~ (89) 

where parameter " a" governs the strength of the mean field attraction and " 6" 
controls a short-range repulsion. Obviously realistic EoS of nuclear matter has 
much more complicated temperature dependence. Due to this we will exploit 
a mVW EoS with 



P[V,T] = f(T)P YW [V,T], (90) 

where f(T) is a function of the temperature. In the given paper we do not 
consider the heat transport, assuming T = const. In this case our solutions are 



self-similar. Doing replacement t — > tJf(T) and (d?7 r + Cr) (<^7r + Cr)/y f(T) 
we recover results for the original VW EoS. 

We choose p T = p cr , T r = T cr as the reference point (see notations in subsect. 
12. ip and perform expansion of the pressure in the vicinity of this point. The 
critical liquid-gas point (evaporation point) is determined from the conditions 
OP/dV = d 2 P/dV 2 = 0, cftP/dV 3 < 0. Critical parameters are 



^~cr — 27^' ^ cr ~ ^Nb, P cr — 271)2' Ucr ~ 

For the given EoS they coincide with those for purely VW EoS. 
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For T = T comp = ||T cr ~ 0.844T cr , n comp = |n cr pressure ( 189]) touches zero. 
For T < T comp there arises density interval (from (| — \J\ — \^) n cr to (| + 
| — |^)n cr ), where pressure becomes negative. In this paper we will restrict 



our consideration by taking T cr . > T > T c 



comp ■ 



For a more convenient analytical treatment of the problem let us, selecting 
corresponding function f(T) additionally fulfill condition §^|(n cr ,T cr ) = 0. It 
allows to parameterize the Helmholtz free energy, as it has been done in the 
paper body for the Landau free energy, with two minima deviating only little 
from each other. Then we may specify function /, e.g., as 



f(6T) ~ C{8T) [1 - 25T] 2 . (92) 

Although our analytical consideration is valid only in the vicinity of the critical 
point < — 5T <C 1 and we perform numerical calculations for < — 5T, let 
us specify the pre-factor C(5T) = [1 + 4(5T) 2 ]~ 1 to reproduce the ideal gas 
EoS for sufficiently low n and high T. 

Expanding the pressure in \8T\ -C 1, 5n/n cr C 1 we get 



QT r)Tr)n QT 

5P = cr / + T7T^(H 3 - QT cr n cr (6T) 2 + ... (93) 
4 Von l CT 

Last term in f[95|) is actually unimportant since addition of any constant (at 
T = const last term is constant) does not change equations of motion. Let 
us count P from its value in the final equilibrium state (reaching at t — > oo). 
Namely this difference P — Pf has the meaning of the thermodynamical force 
driving the system to the final equilibrium state (see (fl3]) ). Then we may 
construct the Landau free energy (fTE]) such that 5(5F L )/S(Sn) = P — Pf. 
Comparing (19"31) with (|T7j) we find relations between coefficients: 



a = 9T cr /(8n cr ), b = l/(3n cr ), v 2 = —3m 2 T cr ST/(2ab) = —A8Tn 2 cr m 2 , 
3ab 9 T cr 

t = 8(d Vv + Cr)(9n cr T cr \5T\y l . (94) 

Since we used f(T — > T cr ) = 1 in Eq. (1921 . these relations are the same as for 
the original VW EoS. 
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7 Appendix B. Mean field and fluctuation region 



In the paper body we considered dynamics of fluctuations assuming that ther- 
modynamical characteristics like pressure, free energy etc., are given (mean 
field approximation) and not modified by fluctuations (mean field approxi- 
mation). In order to estimate a possible influence of fluctuations on thermo- 
dynamical characteristics of the equilibrium uniform system at T 7^ let us 
compare mean field and fluctuation contributions to the density of the specific 
heat 5c v = -T[d 2 5F j {dT) 2 ] v . Using fl22J, (JMD for T near T cr we find 

c^ F ~ 9n cr /2 (95) 
for configurations corresponding |e| <C e max . 

The fluctuation contribution to the specific heat density c' v can be found with 
the help of the functional integration 

exp[5F'/T] = J D5n'exp(5F'[5n'}/T), (96) 
where following ffl6|) we have 

5F' = -j j^y 3 Hk 2 + <y}, a = 2\v 2 /c. (97) 
From here using (|28l) . ( 1941) we obtain 

a 3 / 2 108 ntTl 



167^37/2 71 al\8T\V r 



(9£ 



Equating (|95|) and (|98j) (Ginzburg - Levanyuk criterion) we estimate the 
Ginzburg number 



Fluctuation region is narrow provided Gi 1 and it is broad for Gi £ 1. 
Similar estimate to Gi < 1 follows from the so-called Ginzburg criterion 
^Z 3 |&F MF | > T cr . In the fluctuation region (T cr > T > Tg) fluctuation effects 
may modify ST dependence of coefficients in Eq. ( I93lh 

Substituting in ( 1991) typical values of parameters for the hadron-quark phase 
transition we estimate Gi > 1.4 (lOOMeV • fm- 2 /^o) , i.e. fluctuation region is 
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broad. In case of the NGL phase transition we estimate Gi ~ 10 (T cr /18.6 MeV) 6 
and fluctuation region is also broad. 

One can construct description of the fluctuation region incorporating effect of 
long-range fluctuations directly in the phenomenological expressions for the 
free energy 5F and the pressure 5P. One can do it with the help of the re- 
placements v 2 -> v 2 en = v 2 (T = 0)|5T| 2 / 3 , and A -> A rcn = A(T = 0)|5T| 2 / 3 
(similarly one incorporates fluctuations in description of the superfluid He 4 ). 
Other coefficients in expressions of sub-section 12.21 remain unchanged. With 
these replacements dependence on \8T\ disappears from the Ginzburg crite- 
rion, and both values Cy F and c v do not vanish for 5T — > 0. The Ginzburg 

- Levanyuk criterion then reads as {^■ ncr J i cr ^j <C 1. Effect of not included 
fluctuations is small provided this inequality is fulfilled. 

Finally it is worthwhile to mention that it takes a long time (typically ~ to) 
to develop critical fluctuations. Therefore if the system passes the fluctuation 
region during a time t evoh as it occurs in a course of heavy ion collisions, and 
tevoi < to, critical fluctuations will not have enough time to develop. Thus 
for t < to it is legitimate to use mean field EoS to describe evolution of the 
system. 
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